Exact solution for infinitely strongly interacting Fermi gases in tight waveguides 
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We present an exact analytical solution of the fundamental systems of quasi-one-dimensional 
spin- 1/2 fermions with infinite repulsion for arbitrary confining potential. The eigenfunctions are 
constructed by the combination of Girardeau's hard-core contacting boundary condition and group 
theoretical method which guarantees the obtained states to be simultaneously the eigenstates of 
5* and Sz and fulfill the antisymmetry under odd permutation. We show that the total ground- 
state density profile behaves like the polarized noninteracting fermions, whereas the spin-dependent 
densities display different properties for different spin configurations. We also discuss the splitting 
of the ground states for large but finite repulsion. 
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Introduction. — The experimental progress in manip- 
ulating cold atoms in effective one-dimensional (ID) 
waveguides [l|, Q has stimulated extensive theoretical 
and experimental study of the ID strongly correlated 
atomic systems. Particulary, the experimental realiza- 
tion of Tonks-Girardeau (TG) gases ^, J] has allowed us 
to study the fermionization of Bose gas in the strongly 
interacting limit. More recently an interacting ID Fermi 
gas with tunable interaction strengths has also been ex- 
perimentally realized which offers the opportunity of 
studying the ID Fermi gases even in the TG limit. To 
understand the physical properties of the cold atom in 
the strongly interacting limit, exact solutions and some 
refined methods capable of dealing with strong correla- 
tions are especially important Ji, [7, l8,, JJ] . In the infinitely 
repulsive limit the many-body state of a TG gas can be 
constructed via a Bose-Fermi mapping [9^1 . Despite its 
long history, the generalization to systems including spin 
degree of freedom is a highly non-trivial problem and 
only recently was tackled [l^, [ll[ . Nevertheless the con- 
struction of the exact wave function for the fundamental 
system of indistinguishable spin-1/2 Fermi gas in the TG 
limit is still lacking despite its great importance Dif- 
ferent from the Bose system whose ground state (GS) is 
proved to be a degenerated ferromagnetic state 0, the 
GS of a Fermi system generally falls into the state with 
lowest total spin value S 8]. As we shall clarify later, the 
mixed symmetry of the spin function renders the gener- 
alization of Bose-Fermi mapping to the spin-1/2 Fermi 
system difficult and very challenging. In this work we 
present for the first time an analytically exact solution of 
quasi-lD Fermi gases with infinite repulsion in trapped 
potentials. 

Model. — We consider a quasi- Id system with N spin- 
1/2 fermions tightly confined in an elongated potential 
trap which is described by an effective ID Hamiltonian 
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gid^5{x^- Xj), (1) 

i<3 



with gid being the effective ID interaction strength [14 1. 



For a harmonic pot entia,!, y(xi) = ^mui'^xf. Despite in 



tensive research 15, 16,[l7|, there has been rarely rigor- 
ous results on the interacting spin-1/2 fermion systems 
except the homogenous Yang-Gaudin model [1, Out 
exact solution in the strong coupling limit for the funda- 
mental spin-1/2 Fermi system will provide a firm touch- 
stone for various approximate methods fl^, 1^ 17 1 and 
also deepen our understanding on the few-body system. 

Construction of exact ground-state wavefunction. — In 
general, one can represent the many-body wavefunc- 
tions in the space-spin form as ^' (xi, ai; ■ ■ ■ ; xn, gn) ■ 
Inspired by the seminal work of Girardeau, the effect of 
an infinitely strong interaction can be reduced to a hard- 
core boundary condition 



* (xi,cri; ■ • • ^xn^cfn) |a 



0. 



(2) 



In addition the Pauli exclusion principle enforces the 
boundary condition * (xi, cri; ■ • • ; xn, cfn) \x,=xj,cn=aj = 
for two particles with the same spins. Therefore the 
general contacting boundary condition for a TG Fermi 
gas can be represented as ^ {xt — Xj) — 0, which is ir- 
relevant to the spin configurations. Now it is straight- 
forward to observe that the wave function, which is com- 
posed of Slater determinant of N — + orthonormal 
orbitals (t>i{x), • • • , 4'n{x) occupied by either component 
of Fermions, fulfills the above boundary condition. Ex- 
plicitly, we have 



'iPa{xu...,xn) = iNl)'idet[(f)J{x^)]lIll■ 



.,N 



(3) 



with (j)j{xi) the eigenstate of the single particle Hamil- 
tonian Hi — —-^-^^ + V{xi). So far the spin part of 
wave function is not considered yet. Since H is spin in- 
dependent, H is commutable with the total spin operator 
S — J2i Si , where Si is the spin operator of the ith par- 
ticle. This implies that the system possesses a global 
SU(2) symmetry and the eigenstates of H are simultane- 
ous eigenstates of and Sz- Thus, only the eigenstates 
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with the largest eigenvalue Sz = S are needed to be con- 
sidered and the remaining eigenstates can be calculated 
from them by the lowering operator S^. In addition, 
the total wave function of N indistinguishable fermions 
has to be antisymmetric under transposition of any two 
particles. 

According to the GS corresponds to the fully filled 
state with the lowest N orbital occupied and excited 
states are generated by occupying higher orbitals. Simi- 
lar to the spinor boson case, the GS is highly degenerate 
in the TG limit due to the different spin configurations. 
Among the family of degenerate GSs, the ferromagnetic 
spin state with Sz ^ S ^ N/2 is a product of all spins 
up which is totally symmetric in permutations. The total 
wave function, antisymmetric under transpositions, takes 
a factorized form ^ = ipA (a;i, • ■ • , x^) ■ ■ ■ Xi(^)j 
where Xi(i) denotes the up-spin and x^i'i) the down-spin. 
For the system with fixed up-spin and down-spin parti- 
cles, the ferromagnetic state with Sz = N/2 — m is also 
totally symmetric in permutations and degenerated with 
the polarized state, where m = N^ and n = N^ ^ N — m 
are the numbers of particles with down-spin and up-spin, 
respectively. So far, only the ferromagnetic state is con- 
structed. An important issue here is to discuss how the 
GS degeneracy in the TG limit is split when gid is large 
but finite, or alternatively, to find the GS which could 
be a good approximation of the true wave function when 
the interaction strength is very large but not infinite. Ac- 
cording to Lieb-Mattis theorem [8|, for finite interaction 
strength, the state with lower S has lower GS energy, 
therefore the GS for the system with fixed n and m is 
the state with S = Sz = N/2 — m. Intuitively, the re- 
pulsive interaction term will contribute a positive energy 
to a state with S < N/2, but it does not contribute to a 
ferromagnetic state with all spins oriented in the same di- 
rection, therefore the Lieb-Mattis theorem seems counter 
intuitive. One can understand this problem by noticing 
that a ferromagnetic state with S = N/2 should occupy 
N different orbits due to the Pauli principle, whereas for 
the state with lower S, the particles with opposite spins 
are allowed to occupy overlapping states and thus lower 
the energy. 

The spin function with S < N/2, described by a Young 
diagram [n,m], is not totally symmetric. Nevertheless, 
we can still represent a wave function formally as a prod- 
uct of ipA and "05: where ips denotes a symmetric func- 
tion composed of linear combination of product of sign 
functions and spin functions. Next we shall resort to the 
group theoretical method to construct ipS- 

Before presenting our result, we first introduce some 
notations of group theory Let Ba = {bi,b2, ■ ■ ■ ,brn} 
be a set of m different integers where 1 < 61 < 62 < • ■ • < 
bm ^ N. The n = N — m remaining different integers 
ai, a2, . . ., a-n, satisfying 7^ bj and 1 < oi < 02 < 
. . . < Qn < N are also determined by the set Bq,. There 
are iV!/(m!n!) different sets Ba- bj = n + j when a = 1. 



Corresponding to a set B^, we define a permutation Pa, 

p { 1 2 ... n n+l n + 2 ... N \ 
\ai a2 ■ . ■ an 61 62 ■ ■ ■ J ' 

Remind that Pi is the identical permutation. The 
left coset of a subg roup Sn ® Sm of Sj\[, where 5*^ 
and Sm are respectively the permutation groups of the 
first n objects and the last m objects, is denoted by 
Pa {Sn ® Sjn)- Introduce Qi = nr=i HiLn+i sgn{xi - 
Xj), with the sign function sgn{xi — Xj) = {xi — 
Xj)/\xi — Xj\, and Qa — PaQi- A spin state with 
S ^ Sz ^ N/2 - TO is denoted by = 
3^1"'™' Zq,, which is a basis tensor of the tensor subspace 
J^'^'r of SU{2) of rank N with the highest weight 
where Zi = • • ■ Xl(»^)X2(?^ + 1) ■ • ■ X2(A^), 

Za = PaZi, 3^1"'™! = and = 

(Ei?es„ R) (Et65„ T) {u"U [E ~ (jn + j)]}, where 
E is the identical permutation and (j n+j) is the transpo- 
sition between j and n+j. Our definition for a Young op- 
erator [isl coincides with that in [Toj . but different from 
that in [20|. The spin states with S ^ Sz = N/2 — m 
based on one definition are the linear combinations of 
those on the other. 

Theorem. The totally symmetric wave function con- 
structed by the product of the sign functions Qa and the 
basis tensors Za is 

{W!/(ri!m!) 
E p. {Qi(3^i"'™l^i)} 

Af!/(n!m!) 

= E {yt'^^Qc} Za. (4) 

a=l 

Proof. Since Qi and 3^1"'™' Zi both are invariant in 
left-multiplying by any element of the subgroup 
the action of Pa is proportional to that of the sum 
over all elements in Sn so that ips is invariant in Sn- 
The last formula in ([4]) is obtained by rearrangement, o 

Now it is easy to check that the ground state 

^' = ^Atps (5) 

with tpA and tps given by ([3]) and ([4]) fulfills all the re- 
quirements of symmetry and hard-core boundary condi- 
tion and are simultaneously the eigenstates of S'^ and Sz 
with S = Sz = {N^ — N]/)/2. As a concrete example, for 
TV = 3 with S = 1/2, = ELi(3Qa - l)Za, where 
the identity Eq=i Qq = 1 is used for simplification. We 
note that our constructed exact solution has the same 
spin structure described by the Young diagram [n, to] as 
the true wave function when the strength is very strong 
but not infinite [1]. Therefore, our result is expected to 
interpolate analytically between the finite-repulsion case 
and the limit case with infinite repulsion. 
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FIG. 1: (color online) The GS density distributions of Fermi 
gas in the limit of infinite repulsion with Ni; = 3 and Ni = 2. 



Density distributions. — The spin-dependent reduced 
one-body density matrices are defined as pa-{x,x') = 
C Jl\f^2dx,'^^x,X')P^^'>'^ {x',X'), where X' = 

(X2,--- ,xn), P^J = (1 ± <7i'V2 with a =T (i) cor- 
responding to -I- (— ) and C are normalized constants 
fixed by the conditions / dxpa(x) = No-- Here the spin- 
dependent single particle densities Pa{x) = pa{x,x) are 
the diagonal elements of the corresponding reduced den- 
sity matrices. The total density is defined as p{x) = 
Pi (x) -I- pi {x) . After some algebra, we can prove that the 
total GS density is identical to the density of a polarized 
N-particle free fermion system which takes the following 
simple form p(x) = J2i l<^i(^)Pi where the summation is 
over the lowest N — n + m single particle states. For the 
ID harmonic trap, the orbital are the oscillator eigen- 
states. In contrast with the ferromagnetic ground state 
where Pa{x) = {N„/N)p{x), there is no a simple expres- 
sion of the spin-dependent density for the general state 
corresponding to Young diagram [n,m]. Nevertheless, 
one can calculate Po{x) directly from the exact ground 
state wavefunctions. In Fig.l, we display the GS density 
distributions of a five-particle systems composed of three 
spin-up and two spin-down fermions. Despite the same 
total density distributions, the spin-dependent distribu- 
tions are apparently different from that of the ferromag- 
netic state. Here the total density profile exhibits five 
peaks corresponding to that the wavefunction is com- 
posed of the lowest five orbitals. However, due to the 
spin-dependent term of eq. (U, the spin-dependent den- 
sity profiles are reorganized so that the spin-up part and 
spin-down part avoid overlapping together and show al- 
ternative peak structure to lower the energy. Although 
the exact many-body wave function is constructed, the 
calculation of the density distribution and momentum 
distribution for a large system remains a difficult task 
due to the time consuming to calculate multidimensional 
integrals. Nevertheless, some robust features are found 



FIG. 2: (color online) The GS momentum distributions riaik) 
for the system with N-^ = = 2. Solid line corresponds to 
spin singlet, dotted line to ferromagnetic state and dashed 
line to the case of free fermion for comparison. 

to be not sensitive to the size. For example, the peaks in 
the spin-up and spin-down density distributions appear 
alternatively and the density distributions also show the 
parity symmetry. 

The momentum distribution can be directly obtained 
from the Fourier transformation of the reduced den- 
sity matrices Uaik) = {2tt)~^ J dxdx'e^''^^~^ ^ pcr{x,x'). 
For the ferromagnetic ground state, we have na-{k) = 

iNcr/N)J2iLo^ with (j)i{k) the Fourier transfor- 

mation of the l-th oscillator eigenstate. Similarly, such 
a simple expression does not generally hold true for the 
spin state with the Young diagram [n,m]. As shown in 
Fig. 2, the momentum distribution displays quite differ- 
ent behavior for the spin singlet and ferromagnetic state. 
Comparing with the noninteracting case of gid — where 
n„{k) = X^fcTo"^ l*^/!^)!^ ' the momentum distribution 
develops a wide tail. 

Comparison with system with large but finite 
repulsion. — Apart from infinite repulsion limit, 
there is no analytical solution available for the harmonic 
potential. However, for the small particle system, we can 
apply the exact diagonalization method [2ll| to calculate 
its GS properties and compare with our analytical result 
in the infinite limit. In terms of the fermionic creation 
and destruction operators a\^ and ai^ of the axial 
harmonic oscillator, we get the many-body Hamiltonian 
corresponding to ([T]) 

H = fiiAJx'^{i+]^)d\^di„ + U '^Iijkid\^a^.^dkia^, (6) 

i,(T ijkl 

where = Ix ^T^dx(t)i{x)4)j{x)4>k{x)(t>i{x) are the di- 
mensionless interaction integrals with l^ = ^/hjrrwjx and 
U = gid/lx- The Hamiltonian © can be exactly diago- 
nalized in the truncated basis of eigenstates of the har- 
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FIG. 3; (color online) (a) Comparison of GS density distribu- 
tions obtained by ED for the system with — 2 and A''^ — 1 
with the analytical result of S = 1/2 state, (b) The energy 
versus the interaction for different spin state. The density 
distribution for the ferromagnetic state (inset). 



monic oscillator and then the GS density can be calcu- 
lated numerically. Fig. [3^ shows the 3-particle state cor- 
responding to the Young diagram [2, 1], compared to the 
results obtained from exact diagonalization (ED) of ([6|) 
with parameter U = IBhuJx- The numerical results are in 
perfect agreement with our exact results, which indicates 
that the limit of infinite repulsion is practically reached 
at U = Ibhujx [21, 22 1 . Instead, as displayed in the inset 
of Fig. [31d, the spin-dependent distribution for the de- 
generate ferromagnetic GS state displays quite different 
behaviors from the true GS of the system with large fi- 
nite repulsion. We also show how the GS energy of the 
state with 5 = 1/2 changes versus the increase of inter- 
action strength in Fig. [3|d. In the large repulsion limit, it 
tends to become degenerate with the ferromagnetic state. 
It would be interesting to compare our results with Ref. 
[2^, where the stability of Fermi gases with the pres- 
ence of attractive p-wave interaction has been discussed. 
In the limit case with absence of p-wave interaction, its 
phase diagram is consistent with our conclusion, i.e., the 
GS is an antiferromagnetic state. 

Experimental realizability and detections. — As the Fes- 
hbach resonances in qausi-lD 2-component Fermi sys- 
tem have been observed ^ , in principle the 2-component 
fermionic TG gas can be realized. For the high- 
dimensional system, the physics in the infinitely inter- 
action limit (the unitary limit) as well as the BEG and 
BGS crossover around the unitary limit have been stud- 
ied extensively. The main challenge for realization of 
quasi- ID Fermi system comes from that both the Fermi 
energy and temperature should be much lower than the 



transverse confining energy. The spin-dependent density 
distribution might be detected within the current tech- 
niques which allows for measurement of specie-dependent 
properties [ill . 

Summary. — We have constructed the exact eigen- 
states of the fundamental system of quasi-lD spin-1/2 
fermions with infinite S repulsion by means of group the- 
oretical method. While the infinite repulsion is described 
by a hard-core boundary condition, the group theoretical 
construction guarantees our wave function automatically 
fulfilling the permutation symmetry and being the eigen- 
statcs of S and Sz. The construction scheme and the 
formula for spin densities are valid independent of the 
trapping potential and the particle number. For large 
but finite repulsion we have calculated the ground state 
for a few-particle system numerically by using the exact 
diagonalization method. The numerical result is found 
to be in excellent agreement with our analytical result. 
Our construction of exact eigenstates is valid even when 
a Zeeman term in the Hamiltonian exists because it does 
not change the nature of the states. 
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